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TECHNICAL NOTE 3682 


TIME CORRELATOR FOR PROBLEMS IN AERODYNAMICS 
By George Tolmie Skinner 


SUMMARY 


An instrument of fairly simple design for measuring time correla- 
tion functions of two stationary random electrical signals Is discussed. 

It is intended primarily for use in problems connected with aerodynaml - 
cally produced acoustic fields but has suitable properties for applica- 
tion to a rather wide range of aerodynamic problems involving turbulent 
fields. It has been designed and constructed with a view to economy and 
simplicity of operation and makes extensive use of the general statistical 
properties of the problems for which it is intended. 

A few experimentally determined autocorrelation functions are given 
in order to indicate the degree of accuracy achieved, and the Fourier 
transform of the autocorrelation function of a random input is compared 
with the power spectrum of the same function. 

Some practical points of general interest are discussed. 


INTRODUCTION 


In the description of random fields, the correlation function is an 
important tool. If u(r,t) and v(r,t) are two random variables of a 
stationary, statistically homogeneous field, then all quadratic mean 
values of u and v are expressible in terms of the correlation function 

\|f(r 2 - r 1 ,r) = u( t ) v( r 2 , t + r) (l) 

where the bar denotes an ensemble average, which, for the types of proc- 
esses to be considered here, will be equal to the time average. That 
is, an ergodic property w il l be assumed. 

The space-time correlation function defined in equation (l) can, in 
many cases, be reduced to a space correlation when there exists a mean 
velocity which is large compared with the fluctuations, for example, when 
turbulence in a stream of fluid Is considered and the rate of decay Is 
small. The large mean velocity U permits the interpretation that the 
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measuring instruments are passed rapidly "through, a pattern fixed in time, 
so "that the approximation can he made that 



This, however, is not always possible. 

Recently much interest has been expressed in aerodynamical ly pro- 
duced acoustic fields. Here the above reduction is not possible as the 
acoustic disturbances are not propagated at the mean flow velocity but 
at the speed of sound in the medium. The mean square pressure at a point 
can be expressed in terms of the space-time correlations of the turbulent 
field producin g the noise, as is shown in reference 1. The technique 
of obtaining space correlations simply involves the use of two measuring 
probes at different locations. Introducing the time delay is the diffi- 
cult part of the problem, and for the remainder of the present work only 
the time correlation of two random variables will be considered. 

In addition to the above example, it is possible that in certain cases 
of phenomena connected with spectra having predominantly low-frequency 
contributions it may well be easier to measure autocorrelation than 
spectrum, since it is occasionally difficult to separate resonance peaks 
from continuous spectra with low-frequency electrical filters. Since the 
spectrum can be obtained by Fourier transformation of the autocorrelation 
function, here lies a further possible use of the time correlator. In 
fact, some work has been done in recent years on the study of heart beats 
by nipans of time correlators. In aerodynamics, some buffeting problems 
fall into the category of. very low frequency phenomena. By the Wiener- 
Khinchin relations, if the power spectrum of u(t) is f(cjo) and the 

autocorrelation function is t( T ) = u(t)u(t + f), then 

r»°° 

f(<n) = § / t( T ) cos ^ ^ T 

0 

and 
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The first essential of a time correlator is some means of delaying 
an electrical signal (since the measurable quantities will ordinarily be 
obtained in the form of fluctuating voltages ) . Two - fundamentally dif- 
ferent approaches exist and may be described as follows: 

(a) The signal enters a "box" continuously and emerges continuously 
but delayed in time. This will be designated as the "continuous" system. 

(b) Periodic samples of the signal are taken and stored in a "box" 
from which they are recovered after some predetermined delay time. This 
will be called the "sampling" system. 

The apparatus to be described operates on the latter principle, and 
the reasons for such a choice follow in the next section. ■ Suffice it to 
say here that the signals to be handled will have nearly symmetrical prob- 
ability distributions, and use is made of this in order to evolve as 
simple a device as possible for application to the current aerodynamic 
problems mentioned above. 

This work was conducted at the Guggenheim Aeronautical Laboratory 
of the California Ins titute of Technology under the sponsorship and with 
the financial assistance of the National Advisory Committee for Aero- 
nautics . 


SYMBOLS 


A 

a 

b = l/c 2 

C 

E 

f(<o) 

h(ffi) 

k 

M 

N 


area of output pulse 
velocity of sound 

interstage coupling time, constant of amplifier 
fixed voltage 

power spectrum of an input 

square of response function of an amplifier 

slope of saw-tooth voltage used for sampling 

local Mach number 

capacity of counter 
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k 


n 


n l 

r l ,r 2 

s 

T 

t 

U,U 

u(t) 

v(t) 

X,Y 

5t 

5t* 


counting rate used in storing a sample 
number stored, representing a sample 
vectors; see sketch, with equation ( 2 ) 
number of amplifier stages 

t ime to fill counter at n counts per second 
time 

mean velocity 

stationary random variables 

additive constants used in deriving equation ( 7 b) 
time-interval representation of a sample 
recovered time-interval representation of a sample 
actual input voltage used in deriving equation ( 7 b) 
delay time 


ilr(r 2 - r^T) 


space-time correlation function 


tJt(t) 

flO 

Subscripts: 

m 

max 


time correlation function 
angular frequency 

value of a quantity after passing through an amplifier 

Tna-ylTrmm 


N 

0 


noise 

quiescent value 


T 


turbulence 
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CHOICE OF SYSTEM 


The continuous system usually takes one of two forms. In the first 
of these., the signal is passed through an electrical or an acoustic delay 
line. In the other, the signal is recorded on magnetic tape (or a mag- 
netic drum) and recovered as the tape passes a pickup head, so that a 
delay is introduced equal to the time of transit of the tape between the 
recording and pickup heads. 

The samp ling system may also take two forms. In one, the sample of 
one signal is converted into a number, which is stored in some sort of 
counter, to be recovered later. In the other, the sample is converted 
into a pulse, of length proportional to the value of the sample, which is 
passed through a delay line. 

To choose a system the shortcomings of the elements involved must 
be discussed in the light of the expected applications. In low-speed 
turbulence, frequencies of a few cycles per second axe to be expected, 
while in aerodynamic noise problems, such as those concerning acoustic 
radiation from boundary layers and jets, it is conceivable that frequen- 
cies approaching 100 kilocycles per second may exist. Correlation func- 
tions will then be expected to have significant contributions with max- 
imum delays up to about 0.1 second in the former case, while in the latter 
case the important features may involve maximum delays of considerably 
less than 1 millisecond. 

Clearly it is not likely that one single instrument could cover 
these ranges without becoming a complicated piece of electronic equip- 
ment. However, by suitable combinations of the various techniques it 
is possible that a relatively simple instrument can be constructed and 
with inexpensive modifications can be made to handle any one of the ranges 
at a time. 

To use a tape recorder for delay one must have an elaborate machine 
having many speed ranges and possibly a frequency-modulation system to 
handle very low frequencies. For this reason, and because of the com- 
plicated compensation systems involved as a result of the frequency 
characteristics of magnetic tape and the associated recording and pickup 
heads, the tape recorder was ruled out, although for certain purposes 
it has been found quite satisfactory. 

The continuous delay system us ing an acoustic line was rejected 
because of size and distortion when long time delays and low frequencies 
are to be handled. For delays up to about 1 millisecond, the electrical 
delay line is certainly feasible, but it is felt that its use would be 
more effective in a sampling system, where amplitude distortion can be 
overcome, as will be explained presently. 
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By the use of various very high speed electronic circuits it is 
possible to take a sample of almost any electrical signal. For the range 
of frequencies mentioned above this presents no insurmountable problem, 
but, as will be shown later, for the types of signals anticipated (i.e., 
signals hav in g very nearly symmetrical probability distributions) the 
accuracy with which each sample is measured is of little consequence, 
only the consistency of operation being important. Because of the sta- 
tistical properties of such signals, it is possible, as will be shown, 
to form an accurate time -correlat ion function even with relatively "inac- 
curate" samp ling . 

Of importance is the fact that, once a sample is taken, it is possible 
by means of a pulse length-modulation system to achieve short (less than 
1 millisecond) delays by the simple expedient of an electrical delay line 
and yet by means of electronic counter storage to "hold" a sample as long 
as desired, a week being no particular problem. 

This flexibility as regards available delay ranges makes the sampling 
system attractive. The disadvantage lies in the fact that for long delays 
involving the use of a counter only one sample can be undergoing delay 
at a time; consequently, only a small fraction of the available informa- 
tion is used and time averaging of the results is more difficult. 

A final consideration involves the method of forming the time aver- 
age u(t)v(t + t) . In a continuous system this usually takes the form 
of a pair of squaring circuits (e.g., thermocouples) carrying out the 
operation 

(u + v) 2 - (u - v) 2 = (u 2 + v 2 + 2uv) - (u 2 + v 2 - 2uv) 

= 4uv 

In a sampling system, on the other hand, the delayed sample arrives 
either as a pulse length or as a number (which, as will be discussed 
later, is synonymous with a pulse length). Hence the sample of u(t) 
arrives at time (t + t) in a form suitable for "gating" v(t + t) for 
a short time which is proportional to u(t), thus producing a sort of 
pulse whose amplitude is v(t + r) and whose duration is proportional 
to u(t). The area of this pulse is then directly given by u(t)v(t + r). 
Time averaging provides the mean value. 

It is hoped that the foregoing discussion indicates that the only 
real practical disadvantage of a sampling system for the problems on 
hand (i.e., where signals having nearly symmetrical probability distri- 
butions are to be correlated) arises out of the wastage of information 
when long delays are required, necessitating somewhat longer averaging 
periods. If this does not give rise to undue difficulty, it seems to 
be a small price to pay for the convenience of wide flexibility in delay 
times and ease of multiplication. 
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On this "basis, the sampling system was chosen and designed in prin- 
ciple. Perhaps not surprisingly, it was subsequently found that such an 
instr ume nt had been built at the Massachusetts Institute of Technology, 
(ref. 2). This latter, however, is a very much more general and compli- 
cated sort of instrument, while the one described herein relies heavily 
on the statistical properties of the expected random fields to be dealt 
with for its accuracy and simplicity. 


PRINCIPLES OF OPERATION 


The overall function of the instrument is to perform the following 
operation on two stationary random variables u(t) and v(t): 

i}t(t) = u(t)v(t + t) 

For the sake of clarity, an approximate treatment of the principles of 
operation will be given first, followed by a more exact analysis. 

Approximate Treatment 

Consider a slowly varying function of time u(t). At time t-^ let 

there be added to u(t) a linearly increasing voltage k(t - tj). If 
the variation of u(t) during this addition is neglected, the sum will 
reach a predetermined level E, sifter a time interval 6t, given by 

u(ti) + k6t = E 

whence 

5t = | “ ^(^1) ^ 




8 


NACA ET 3682 


If a pulse generator producing pulses at the rate n per second, is 
connected, to an electronic counter of capacity If by a "switch" which 
closes only during St, the operation may he denoted hy 



n/sec St Up 


where 


np = nSt 


(5) 


Clearly, St does not have to he an integral number of times the 
pulse spacing l/n, hut if n is assumed large enough equation (5) is 
very nearly exact. 

The delay time r may he generated hy picking off a point on any 
monotonically rising voltage which starts at time tp. For example, a 
slow saw tooth can he used. The time r can he varied hy raising and 
lowering the pickoff voltage, and at time (tp + r) a signal, such as 

a sharp pulse, can he generated to i ni tiate the multiplication of u(tp) 
hy v(tp + t). 


t-i tp + t 

b* r H 



This process is carried out as follows: 

At time (tp + t) the switch between the pulse generator and counter 
is again closed, so that digits are fed into the counter at the rate n 
per second until the counter fills up to its maximum capacity H, at which 
time it resets itself to zero and opens the switch, turning off the supply 
of pulses. The time taken for this filling-up process is given hy 

St* - - Up) 

For convenience, a t ime T may he defined to he the time it would 
take to fill the counter, starting from zero, if the pulse generator were 
continuously connected to it. Clearly, T = N/n, whence 


6t* = T - 5t 


( 6 ) 
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Suppose, now, that there is an output terminal which is connected 
to v(t) only during St*. Neglecting the variation of v(t) during 
St*, this terminal will carry a pulse of approximately rectangular shape 
whose amplitude is v(t^ + r) and whose duration is St*. The area of 
this pulse, then, is given by 

A = $t*v(tq + t) 


Using equations ( 6 ) and (4), 


A = 




u(ti) 


vit], + 


t) 



Finally, if A is averaged over many repetitions of the process 
and it is assumed that v(t) =0, then 

A = 5 u(t)v(t + t) = ^ t|t(t) (7) 

If u(t)=v(t), then the autocorrelation is obtained. 


More Exact Treatment 

Having established the principles of operation, a more detailed 
anal ysis will be followed. Here the essential difference will be that 
the variation of u(t) and v(t) will be tahen into account by including 
the linear terms of the Taylor expansions of u and v about the points 
t = tq and t = t]_ + t, respectively. Making use of the fact that the 

distribution functions of u and v are nearly symmetrical, it w ill be 
shown that the first correction term in equation ( 7 ) involves the second 
derivative of \Jr( t) . 

With the same numbering system for the equations as was used in the 
section "Approximate Treatment," 

u(tq + 5t) + kSt = E 
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whence, writing u(tp + 6t) = u(tp) + 5tu'(tp) and assuming k » u'(tp) 
(this is a practical requirement). 



u(tp) 

k 


“ ^ u '( t l) 


. ^(tpju^tp) 

k2 


(M 


Let 


n x = n(St - e) (5a) 

where 0 - 6 < ^ and where n^_ is an integer. The number registered 
during 8t may then he np with prohahility 1 - ne and np + 1 with 

prohahility ne. That is, there is an uncertainty of one count in the 
number stored during 5t. This is akin to the problem of measuring the 
distance between two points with a ruler having very few graduations. 

If, for example, the points are 2.7 inches apart and the ruler is grad- 
uated in inches, the number of graduations appearing between the points 
will be either two or three. Now it is easy to show that, if the ruler 
is laid across the points many times at random, the numbers two and three 
will occur with such probabilities that the average result will be 2 . 7 . 

In fact, if the graduation marks have a width w, and a graduation is 
counted if any part of it appears between the points, then the average 
result will be 2.7 + w. The width w corresponds to finite width 
pulses and merely adds a constant to the result. 


Since the counting rate n is very high compared with the rate of 
sampling and these two are not in any way synchronized, it may be said 
that the interval 5t is laid on the counting pulses at random. If only 
linear expansions of u and v are to be considered, it may then be 
said that the expected value of np is n 5t and the expected value 

of 5t* is given by 

5t* = T - 8t (6a) 

The area under the v(t) curve presented during 5t* will then be 

p tp+T+St* 

A = / v(t)dt 

' J tp+T 


and, using the linear expansion for v(t), 

p tp+T+5t* j- -1 

A. - / v(t x + t) + (t - tp - T ) V * ( tp + t) 

■ J tp+T L J 


dt 


= v(tp + -r)8t* + v'(tp + t) 


(5t*) 2 

2 
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Hence , 


A 


v(tp + r)8t* + 


v'Ct! + t) 


(8t*) 2 

2 


(7a) 


which is evaluated in the appendix, for the types of functions to he 
considered, giving the result 


A 


g jConstant + ^(t) + 0 


(7b) 


The output circuit of the correlator is actually a balanced circuit 
in which the constant term is eliminated. 


As a test, the correlator was checked on the autocorrelation of a 
square wave, since this has discontinuities in \jr(r). This will be dis- 
cussed later. 


DESCRIPTION OF CIRCUIT 


The circuit consists of chains or loops of relatively simple ele- 
ments, each of which performs one of the operations described in the 
section entitled "Principles of Operation." The instrument incorporates 
slightly more than 40 vacuum tubes, mostly 6 AK 5 miniature pentodes, 

12AX7, 12AIJ7, and 5963 twin triodes, and 6AL5 twin diodes. These are 
employed in such well-known circuits that little will be said about them ' 
and no circuit diagram will be attempted. The only critical adjustments 
are readily made with the aid of a good direct-coupled oscilloscope (which 
incidentally, is part of the equipment, being used for setting the delays) 
As the critical adjustments depend on the exact mode of constriction, a 
circuit of the present setup would be useless for any other arrangement. 
Considerable thought went into the arrangement of the circuit loops in 
order to insure the fastest possible operation of the elements, and, for 
example, the diode-coupled 59^3 scale-of-two counter storage operates 
very satisfactorily at 300 kilocycles per second. 

A point perhaps worth mentioning is that, in this circuit, the 
flip-flops are driven through a separate diode for each operation; hence, 
in the first place, they cannot become confused and, further, their 
transitions follow the sharp edges of the triggering pulses. The fast 
rises at the plates are never restrained by the tail of a pulse. Twin- 
triode flip-flops were found perfectly satisfactory. Voltage discrim- 
inators, on the other hand, were built from pairs of 6AK5 vacuum tubes. 

To avoid confusion,, the circuit elements will be defined, and the 
symbols used in the schematic diagram (fig. 1 ) are as follows: 
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(1) ADD, adding, network: Passive-resistance network giving an out- 

put proportional to tie sum. of two or more inputs. 

(2) AV, averager (or integrator): A device which forms the time 

average (or integral) of the input signal. 

(3) C, counter: Chain of scales-of-two (six in this case) which 

puts out a pulse when it is full (64 counts). 

(4) D, discriminator: Circuit having two stable states, determined 

•by the input voltage being above or below a predetermined value. The 
circuit puts out a pulse at transition. 

(5) FF, flip-flop: Circuit .having two stable states which is trig- 

gered by an input pulse. Output voltage has two possible values. As 
used in this setup, the inputs are separated. 

(6) GA, gated amplifier: Amplifier to which v(t) is applied and 

which gives an output of zero except when it is gated open during St*, 
when its output is proportional to. v(t). In the circuit it is actu ally 
a difference amplifier having inputs v(t) and -v(t) and is gated via 
the common cathode resistor. 

(7) IB, time base (or linear sweep generator): Circuit which pro- 

duces a voltage increasing from zero linearly with time. It Is started 
by an input pulse and stopped after a time determined by a gating tube 
within the circuit. 

(8) V, dire ct - current voltage source: Potentiometer to give adjus- 

table steady voltage for setting signal levels. 

(9) M, mixer (or coincidence circuit): Circuit producing a fixed 

output when both inputs coincide — acts as a switch between the contin- 
uous pulse generator and the counter. 

(10) 0, oscillator (or pulse generator): 300-kilocycle-per-second 

oscillator followed by an amplifier which suppresses the negative halves 
of the waves. This is the counting frequency n per second. 

The operation of the instrument centers around the flip-flop FF]_. 
(See fig 1. ) An oscillator supplies the sampl i ng frequency, and a dis- 
criminator D]_ picks off a point on each cycle and from it produces a 
pulse which initiates one cycle of operation. Both time bases are trig- 
gered by this pulse, and the flip-flop EF]_ is tripped, causing the 
mixer M to feed the counting frequency from the osci llat or 0 to the 
counter C . 
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The fast t ime -base output is added to the input u(t) and when the 
Siam reaches a fixed value E the discriminator Dg retriggers FFp, 

cutting off the oscillator 0 from the counter C and leaving C with 
a certain number registered. 

Meanwhile, the slow time-base output continues to rise, and when 
it reaches a preset value, determined by the triggering level of the 
discriminator D* and the additive voltage V, the delay pulse is gen- 
erated, once again triggering FF]_ and supplying the counting frequency 

through M to C. Flip-flop FF 2 follows FF]_ through this part of 
the cycle. 

This state continues until the counter is "full," and, as it returns 
to zero, it sends out a pulse to return both FFp and FF 2 to their 
initial states. 

So far, FF-j_ has gone through two cycles, and it will be observed 
that FF 2 follows it only on the second cycle. Flip-flop FF 2 drives 

a cathode follower whose load resistance is also the common cathode 
resistance of the difference anplif ier GA. Ordinarily GA is cut off, 
but when FF 2 follows FF-j_ on its second cycle, GA allows v(t) to 

pass, thus producing an output approximately described by St*v(t + t). 

Finally, averaging, or integrating over a period of time, forms 
\|r(r) as described in the theory. 

Averaging was first accomplished by connecting a long-period ballis- 
tic galvanometer across the plates of the difference amplifier GA. The 
plate loads of the amplifier were bypassed with large capacitors in order 
to effect partial averaging before taking the output to the galvanometer. 
This arrangement proved quite satisfactory for periodic inputs but was 
wholly unsatisfactory when tried on autocorrelation of a turbulence signal. 

The next step taken was to construct GA from a pair of 6AK5 pen- 
todes operating on 90 volts, with the load resistors remioved and purely 
capacitive loads Installed. A bias of 20 volts would reduce the plate 
current to about l/lOO microampere so that little leakage took place 
during "off" periods. By charging the plate-load capacitors repeatedly 
for 30 seconds and discharging each time into a ballistic galvanometer, 
it was found that good averages could be obtained with from 5 to 10 min- 
utes total integration time. In order to do this, 10-megohm plate resis- 
tors were used with 0.1-microfarad capacitors bypassing them to a point 
held at the mean plate voltage (to reduce leakage variations).. A KA-O-ohm 
resistor and 0.1-microfarad polystyrene film capacitor were connected in 
series across the plates, and the voltage across the capacitor was taken 
to an electrometer tube operating as a cathode follower and driving a 
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microarmnet er . This latter part is essentially the "microvolter" described 
by the Victoreen Instrument Co. in its electrometer-tube catalog. The 
microvolter and the battery for the galvanometer lamp are set on paraffin 
blocks and al 1 connections are mode with cotton-covered wire boiled in 
paraffin wax. 

The integrat ing resistor and capacitor are actually incorporated 
in the microvolter alongside the electrometer tube, The time constant 

of the system is 15 ^ minutes and either 5- or 10-minute integrations 

are used. The results presented here for random input were taken with 
5-minute integrations. If the same time is used for each reading, the 
results can be used directly, since some form of normalization will 
determine the scale. 

The operation of the integrator was so consistent that it was used 
also in the spectrum measurements taken to check the autocorrelation 
curve discussed later. The output meter of a Hewlett-Packard heterodyne 
wave analy zer was replaced by a 2,200-ohm resistor (to increase the out- 
put voltage). The integrator was connected across the resistor in 
series with two 5-megohm resistors and allowed to integrate the voltage 
for 5-minute periods, completely eli min ating the dubious procedure of 
trying to guess the average position of the dancing meter needle. 

Possible modifications to the correlator as a whole might involve 
the replacing of the counter storage by a delay line for very high fre- 
quency signals, at which time the sampling process can be speeded up. 

More accurate representation of the sample would be possible without the 
digital storage and more than one sample could be present in the line at 
one time. 

For very short delays, with the instrument as it stands (e.g., for 
the study of correlations near r = 0), a system of interlaced counting 
pulses derived from a square wave could be used to allow St* to overlap 
St to some extent and hence remove the requirement that St be com- 
pleted before St* can start. 

For very long delays, each individual product u(t^)v(tj[ + t) could 
be integrated and the results accumulated in a digital memory. 

Several variations are undoubtedly possible to cover various require- 
ments as they arise. In each case, of course, the limitations of the 
instrument must be considered. 
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SOME PRACTICAL EETAHS 


In the section entitled "Principles of Operation" a quantity A 
was defined. It is the output of the correlator per cycle of operations 
and may he thought of as an essentially rectangular pulse of .current. 

There are certain physical limitations to its size. Its average duration 
will he half the fill-up time T of the counter. With a counting rate 
of 300 kilocycles per second and a storage capacity of 64, T/2 is of 

the order of 10 - ^ second. If the rate of sampling is about 100 per second, 

then the "duty ratio" will be about 10“ 2 . That is, the overall average 
current will be about l/lOO of the average current per pulse. To give 
some values without going into the details of the choice, the output 
pulse is drawn by a tube from a load consisting of 10 megohms bypassed by 
0.1 microfarad. For a mean voltage drop of 100 volts across this resistor, 

the mean current will be 10~5 ampere. The average current per pulse, then, 

will be 10“3 anfpere or 1 mill! ampere, which is a reasonable value consid- 
ering that the tube must be operated at fairly low voltages in order to 
obtain small leakage when it is cut off. 

By using a balanced-circuit arrangement this mean drop of 100 volts 
is canceled out, and it is of interest to inquire by how much it changes 
when the correlation coefficient goes from 0 to 1. The situation may be 
visualized as shown in the accompanying sketch. 
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The solid rectangle shows the maxi mum size of a pulse, its minim u m size 
being zero. With zero fluctuating input, the shaded rectangle Aq is 
obtained. The distributions of the signals u and v are drawn on the 
time and current axes, respectively, and the me anin g of the a ddi tive con- 
stants mentioned in the theory should now be apparent. Since only the 
fractional chang e in Aq is sought, no loss of generality arises from 

choosing eq ual scales for time and current. 

Suppose that v = u and no time delay is introduced. That is, 

the instr ume nt is to form \jr(o) = u(t)u(t) = u2(t ) (which, incidentally, 
it c anno t do, since it is necessary that t be greater than Stmax)* 

This would represent the maximum change in output expected. 

Now in turbulence it is a fairly general result that the peats of 
the signal lie within four times the root-mean-square value. Hence the 
gain controls on the inputs can be adjusted so that the peaks just cause 
the rectangle to run between the maximum value and zero. Hence each 
rectangle is given by 



when v = u and r = 0. Hence A = Aq + u2 and, from the above state- 
ment, \/u2 = h/8 so that u2 = h 2/64 while Aq = h 2 /4. Consequently, 

A-Ao_u 2 _4_1 

“ 25 “ 65 “IT- 


For linear operation it is necessary to have some extra current 
per pulse -to insure that v(t) never causes the tube to operate too 
close to cutoff during the "on" period. So, in practice, although the 
above estimates would indicate a mean change in plate voltage of 
100/l6 volts, it is more realistic to expect only about 3 volts. For 
a balanced circuit the change is doubled, with the result that the unity 
correlation coefficient of a turbulence signal corresponds to a mean 
output of something like 6 volts. This is sufficient to give good accu- 
racy with the integrating circuit. 


A USEFUL TRICK 


In the estimate carried out in the section "Some Practical Details" it 
was ass ume d that the peak of a turbulence signal does not exceed four times 
its root -mean- square value and it was assumed that the signal level can be 
adjusted to the point where such a peak value will just fill the counter 
storage. This is only appro ximat e, and an occasional peak will cause the 
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counter to overfill, with the result that at a count of 64 it returns to 
zero and starts again. Consequently, 6t*, which should he negative in 
such a case (an impossibility, of course), comes out rather large and 
positive. If many peaks which give rise to this phenomenon occur, con- 
siderable error will result and it would be more accurate to clip the 
signal so that such a high peak would lead to 6t* = 0. 

In reference 3 it is shown that for cross correlations of Gaussian 
signals, if one signal is clipped, no matter how severely, the correlation 
remains correct except for a constant factor. Turbulence distributions 
do not deviate far from Gaussian distributions; hence, a certain amount 
of peak clipping should be tolerable. Such a step leads to larger outputs 
from the instrument for the same correlation function, and, up to a point, 
the accuracy is improved. 

Ahead of the u(t) attenuator, a box containing four 1-i volt cells , 

two germanium diodes, and suitable resistances was ins erted in the line 
to limit the peaks to approximately +3 volts. The amount of clipping is 
determined by the gain of the amplifier supplying the signal, and an 
attenuator in the correlator adjusts the result to suit the counter storage. 
The check result on turbulence was carried out using this crude clipper. 


TIME CORRELATOR IN PRACTICE 
Effect of Low-Frequency Cutoff 


In the present section the relationship between the autocorrelation 
of a fluctuating velocity component in a turbulent fluid flow (e.g., 
downstream of a grid) and that obtained from the correlator by analyzing 
the signal delivered by a properly compensated alternating- current hot- 
wire amplifier will be discussed. The following quantities will enter 
the discussion: 

The turbulence spectrum f(a>) will be such that 



f(u>)ckD 



and f(o) i 0. 

The autocorrelation of the fluctuating velocity u will be 


i}t(t) = u(t)u(t + t) 
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These are related by 


r “ 

f(co) = / t[c( t ) cos air dr 

and 


= / f ( cd)cos odt da 

J 0 

It may be noted that iJf(o) = u.2. 

Let it be assumed that the time correlator is perfect. It ■ ft 11 
then form the exact autocorrelation of the signal it receive, from ;he 
hot-wire amplifier. This signal -will, in general, differ from the tur- 
bulent velocity in some important respects, primarily as a result of the 
low-frequency cutoff of the hot-wire amplifier. 

If the output U£ of the amplifier is related to its input u-^ by 
the relation 

( S |) 2 - *-> 

where h(o) =0, then the effect of the amplifier may be expressed in 
terms of h(<n) . 

Using subscript m to denote "measured" (or "entering the corre- 
lator"), 

f m (a>) = h(cD)f(co) 

so that 


Tjr m (-r) = / h(Go)f (cd)cos ojt du> 

m Jo 

and 


% 2 = t m (0) = P h(©)f(a>) dn 

Jo 

The amplifiers common in turbulence work have low-frequency responses 
which are level down to less than 1 cycle per second; consequently, there 

is not much difference between 1 %^ and u^, or between i(r m (o) and ijr(o) . 
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However, since h(0) = 0, so also f m (o) = 0, with the unf ort unat e 
result that, while 


2 

2 


■ 00 

f(0) = P ^(t) dr 0 
J 0 


for the actual turbulence. 


f f m (o) = f V( T ) dT = 0 


for the measured signal. 

There is, of course, some argument as to whether f(0) = 0 in a 
wind tunnel of finite size, but certainly the cutoff occurs at fractions 
of 1 cycle per second much smaller than the amplifier cutoff. 

Assu min g a power spectrum which approaches f(o) horizontally, 
the effect of an amplifier whose low-frequency cutoff is essentially 
contained within the region where f(oi) « f(o) can he derived as follows 


*(0) - i|r m (°) = 



1 - h(co)Jf(tt)) 


cko « f(o) 



1 - h(co) 


d£D 


It follows that 


»(0) - * m (0) 
f(0) 


is approximately equal to the fractional 


loss of energy due to the low-frequency cutoff of the amplifier. This 
is usually negligible in hot-wire anemametry equipment. 


To obtain the effects for large values of t the computation can 
be carried through fairly easily for the case of a "white" spectrum, 
and this shows the essential features. The amplifier may be taken to 
be one having s stages, each involving a resistance-capacity coupling 
of time constant C. That is, the coupling time constant is t he same 
for all stages throughout the amplifier. Then 


f(co) = f(0) = Constant 
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t(*r) - tm( T ) 

~f(o) 


■ X" f - frffe)' 
■ri‘-(A)' 


COS CUT dCD 


COS CUT < 3 cu 


where b = l/C 2 . 

Integrating by parts: 


1 

£ 

1 

Ti 1 ^ f\ 

sin cut 

2b r a^- 1 

f(0) 

L U + ^/J 

T 

0 T 0 (b + cu^) S+1 


sin cur dcu 


Now, as cu — >°°, 


CD 


b + CD 2 

the first term vanishes at both limits. 


1, and, as cu — >0, sin cur — >0 so that 


Therefore 


* - *m = a f “ 

f(0) T ^0 (h +(D 2) s+1 


sin cut dcu 


This Fourier sine transform can be found in reference 4. It is 
necessary that 0 S (2s - 2) = 2s, which is automatically satisfied, and 
that |arg b | < it. Since b = l/T 2 is real and positive, the latter 
is satisfied also. 


Then 


*_-Jk = 2b- (-I) 28 " 1 * Jl£_ ( b s-l e -bl/2r) 

f( 0 ) T s« 2 dljS \ ' 


= -Jtb d^_ ^s-le-bl/Sr) 
TSI ,ThS \ * 


dbs 

Using the expression for the sth derivative of a product, this becomes 7 


si 


- ^m _ J^tb y _ 
f(0) rsi (s - k)i kl 


db 


= ZH* y b k ~l (s - 1)1 d k ( b 1 / 2 ^ 

t ^ (s - k): ki (k - 1 ): abh \ > 
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which is of the form 

-^(0^ = ^ e" bl ^ T (o^r + c^r 2 + . . . + a s r s ) 

= -jtb (oj_ + 027 + . . . + o s t s “^ 

By substituting x = br^ it can be seen that the number of zeros of 
this function is the same as that of f*), which is (s - l). 

Hence - will cross the t axis (s - l) times and will 

approach zero as r— This feature is noticeable in some of the 
results in references 5 and 6. 


For the simple case of a single dominating time constant, s = 1 

and 

V 7 fa - “flb d / bl/^r) _ * • b l/2__ b l/ 2 T _ it p-r/c 

f(0) T dFv ' ~ 2 6 “2C e 

which is a monotonic decreas in g function of t. It must be remembered 
that this was computed for a white spectrum and hence is not valid 
near r = 0. 


Separation of Acoustic Radiation From Turbulence 

In a turbulent fluid motion that part of the fluctuation field 
which travels at the mean stream velocity may conveniently be defined 
as turbulence, while the disturbances which travel at the local sound 
velocity may be considered as noise. When this is done, the two parts 
can be separated, assuming that the acoustic field consists of plane 
waves traveling in one direction, as follows: 

Consider two hot-wires placed a distance Ax apart and alined 
with the mean flow velocity ’ U. Let there be plane acoustic waves 
traveling at velocity a also in the direction of U, and assume that 
the rate of attenuation of these waves is small. 



The fluctuation velocity consists of two parts: 
u(x,t) = Urji(x,t) + Ujj(x, t ) 
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where the subscript T denotes turbulence and N denotes noise. Then 

u(xi,t) = uj^^t + + Ujj^Xgjt + 

where the decay of Urp is assumed negligible over Ax. 

Let 

i1t(t) = u^jtju^jt + r) 


= U T (X2,t)U 3 ,^X 2 ,t + T + + u N (X2,t)u N (x 2 ,t + T + + 


Urp^t^UN^t + T + + U N (x 2 ,t^U T ^X2,t + T + ^ 

= *tt( t + if ) + %n( t + ¥) + %n( t + ¥) + %r( T + Jr ) 

If u T and ujj are independent, then TjfrjTjj = = 0. If Ax is large 

enough so that then ijr^ t + 0 ^ or 811 values °f 

t — 0, so that, approximately, 

'Kt) = %r( T + Tf) 

Alternatively, it may be said that, if the autocorrelation of Up 
is zero for delays greater than or equal to Tq, Ax must be chosen 
such that Ax = Utq. Then i(f( T ) = tyjjpj( T + Mtq ) f where M is the local 
Mach number. 

As an example, consider isotropic turbulence at atmospheric pressure 
with free-stream U = 1,000 centimeters per second. Producing the tur- 
bulence with, say, a l/2-inch grid, it would have a characteristic time 
of about l/2 millisecond. So at a 10-millisecond delay the autocorrela- 
tion should be sensibly zero. That is. Ax should be about 10 centimeters. 

Since a = 33 X 10^ centimeters per second, M = 0.03 and TqM = 0.3 
millisecond. Delaying either u(x^_,t) or u(x 2 ,t) enables the whole 
autocorrelation curve of the noise to be obtained. 
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PERFOBMAHCE 


The correlator, as it stands., has heen tested on sine waves, square 
waves (which, of course, give discontinuities in the derivative of \}t(t)), 
and on turbulence behind a grid. These results are shown in figures 2, 3 # 
and Since the correlator is not designed to take square waves, the 
discrepancies near the negative peaks in figure 3 are most likely due to 
asymmetrical distortion of the wave form. 

Figure 4 shows the result after a turbulence signal has been passed 
through an amplifier having one dominant short time constant. The exper- 
iment was actually carried to much larger values of t, and, within 
experimental accuracy, a second .crossing of the axis could not be detected. 
It thus bears out the features discussed in the last section. 

Figure 5 shows the Fourier transform of figure plotted along with 
the measured power spectrum. 


CONCLUDING REMARKS 


When the effects of auxiliary equipment are not forgotten, the time 
correlator appears to be satisfactory for work in turbulent flows. For 
boundary-layer studies of the type which might arise in acoustic radiation 
problems some modifications may be necessary. By defining a time scale 
from a characteristic length and the flow velocity, an idea of the signifi- 
cant time delays can be obtained. The grid size in the plot of autocorela- 
tion of turbulence was about l/2 centimeter and the flow velocity, about 
10 meters per secondj hence, the time scale may be taken as l/2 millisec- 
ond. The maximum delay required is about 10 times this value. Carrying 
out the same estimate for a boundary layer on a flat plate 10 centimeters 

long at a Reynolds number of 5 X 10^ and forming a time scale from the 
boundary-layer thickness and the flow velocity, one obtains a value of 
about 100 microseconds. Maximum delays of about 1 mi lli second might be 
required. If so, the rate of sampling could be increased five times, while 
the speed with which each sample is taken could be increased correspond- 
ingly, so that the details of higher frequency components could be caught. 

That is to say, as the phenomenon contains higher and higher fre- 
quencies, the time scales should go down, and the appropriate c h a n ges 
to the correlator fit it to the task. 

Instrumentation in the way of probes to pick up such phenomena is 
under continual development. 


California Institute of Technology, 

Pasadena, California, July 11, 1955* 
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' APPENDIX 

DERIVATION OF EQUATION (7b) 

The derivation of equation (7b) is given below: 

A '= v(t + r)8t* + v'(t + (7a) 

Substituting from equation (^a) into equation (6a), 

8t* = (T - E/k) + (l/k)u(t) + (E/k 2 )u'(t) - (l/k 2 )u(t)u' (t) 
or 

k8t* = (kT - E) + u(t) + E U '^ - u(t) ^ 

Squaring and neglecting terms containing 

k^Bt*) 2 = i (kT - E > 2 + i u 2 (t) + (kT - E)u(t) + (kT - E)E - 

(H - E)u(t) + Eu(t) - u 2(t) 

Henceforth it will he understood that u and u' are evaluated 
at t, while v and v* are evaluated at (t + r). So, 

kA = (kT - E)v + uv + E ^ - 32^ + 

i(ET-E) 2 H 1 Si «. (KT - E) 

where second-order terms in and have been omitted. 
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Since u and v axe stationary, random signals, 
tIt-jjlCt) = u(t)v(t + t) = u(t - r)v(t) 
therefore 

^'(r) = u(t)v'(t + t) = -u'(t - r)v(t) = -u'(t)v(t + t) 
and 

i|r 21 (r) = u 2 (t)v(t + r) = u 2 (t - T)v<t) 


therefore 


i|r 2 l ’( t) = u 2 (t)v*(t + t) = -?2u(t - t ) n 1 ( t - r)v(t) = -2u(t)u'(t)v(t + r) 
That is, uv* = -u'v, u 2 v' = -2uu'*v, and v’ = 0; hence 


kA = (kT - E)v + uv 4- (2E - kT) ^ - 2 3^^ 

From the principles of operation it should he clear that u and 
v are always positive. The actual inputs, however, have zero mean 
values and sufficiently large constants are added to these in the instru- 
ment. let the actual input voltages he £ and tj, such that | = ij = 0, 
and let the additive constants he X and Y. Then u = X + £, 
v = V + T), and 

kA = (kT - E)Y + XY + gi] + (2E - KT) ^ - 2X 4^ - 2 AO 

K. K. 


= Constant + §Ti + (2E-kT - 2X) A_H _ 2 Al-Jl 

jS. j£ 

In operation, the adjustments are such that, approximately, kT = E 
and 2X = E. 

So (2E - ME - 2X) « 0 and the presence of g’/k in the correla- 
tion makes this term of smaller order than Jr]. In reference 7 it is 
shown that for autocorrelation (g = tj) this term would he identical l y 
zero in the case of turbulence. Reference 8 shows that, for random 
variables having symmetrical probability distributions, odd-order cor- 
relations are zero. Turbulence, of course, is not quite symmetrical in 
the derivatives of the fluctuating velocities, but the deviation is not 
large. Hence the last term can be assumed negligible, particularly as 
k appears in the denominator here also. 
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Hence, to first order, one obtains 


A = 


1 

k 


[c 


Constant + 



■where -^(r) = g(t)i](t + t). 


Now that the additive constants have been absorbed into one constant, 
u and v may be considered as the actual inputs and the notation 


Tjr(-r) = u(t)v(t + t) 

may be resumed where it is now understood that u = v = 0. 


It is not hard to see that a second-order correction will give rise 
to terms of the form uv"/k2, u'v'/k^, u"v/k2 , and so on. 

These are the second derivatives of ^(t), so that the expression 
for A may be written 


A 


i •! Constant + t ) 


l 



or, since 


E/k « 


A = i j Constant + t( T ) + 



(Tb) 
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Figure 1.- Block diagram. 










Figure 3 .- Autocorrelation of a square vave. 
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Figure 5.- Check against spectrum. 
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